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Abstract 

We study a Hele-Shaw problem with a mushy region obtained as a 
mesa type limit of one phase Stefan problems in exterior domains. We 
study the convergence, determine some of the qualitative properties and 
regularity of the unique limiting solution, and prove regularity of the free 
boundary of this limit under very general conditions on the initial data. 
Indeed, our results handle changes in topology and multiple injection 
slots. 
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1 Introduction 

Given a bounded domain D 6 IR n with smooth boundary dD, a finite set 
of closed curves {s J (x) = 0} such that D is contained in the union of their 
interiors, and a continuous function p(x, t) defined on dD x (0, oo), the classical 
Hele-Shaw problem with Dirichlet data is frequently formulated as follows: 
Find a function V(x, t) and a family of domains S(t) (which each contain D) 
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with dS(t) = {t = s(x)} and such that 
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(l.l) 



This problem models the advance of the slick formed by injecting oil be- 
tween two nearby plates, and has further been used in injection molding (used 
in turn in the packaging industry, and more generally for the production of 
plastic components, for example interior pieces of cars and aircraft), in elec- 
trochemical machining (see |MR] ) . and even to predict tumor growth (see 
[BE])- Sometimes the normal derivative of V at the "slot," dD, is prescribed, 
or curvature dependent terms are included in the free boundary condition. 
Among the most pressing open questions concerning the Hele-Shaw problem 
are finding a weak formulation, studying the regularity of V(x, t) in t, and 
determining the regularity of the free boundary t = s(x). Another question 
which has long attracted interest is whether a Hele-Shaw problem could have 
a "mushy" region. 

In turn, the "Mesa" problem describes the limit pattern lim m ^oo u m of 
solutions u m of, say, the porous medium equation, when the initial data are 
held fixed. This problem first appeared in connection with the modeling of 
problems related to transistors (see [EHKO ). Caffarelli and Friedman studied 
the initial value problem in H n x [0, T] in |CF| . They proved that the limit 
exists, that it is independent of the chosen subsequence, that it is independent 
of time, that it is equal to the characteristic function of one set plus the 
initial data times the characteristic function of the complement of that set, 
and finally that that set can be characterized as the noncoincidence set of a 
variational inequality. Further developments showed that the same conclusions 
hold for the limit when u m is replaced by a fairly general monotone constitutive 
function <f>{u) with 0(0) = (see |FHj for example), moreover, this behavior is 
a property of fairly general semigroups (see jlj, |BEG| ). A Mesa problem for 
an equation which gives rise to a mushy region was studied in |BKMj . For the 
Mesa problem we study in this paper we show all of the properties shown by 
|CFj mentioned above, except that our limits will not be independent of time. 
(This evolution in time is natural since we work in an outer domain where the 
inner boundary data serves as a source.) 

In this paper the authors exploit a Mesa limit setting in an outer domain 
D c to obtain naturally a weak formulation of the Hele-Shaw problem (with 
Dirichlet condition as in Equation ( 11.1)) on the slot as above). The use of 
one-phase Stefan problems with with "mushy" regions and with increasing 
diffusivities naturally produces a mushy region when we permit initial data 
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m for the approximating problems to take values in the interval [0, 1]. These 
Ui can be thought of as generalized characteristic functions. Another aspect 
of this approach which is extremely attractive is the fact that changes in the 
topology of the "wet" region do not interfere with the construction. In short, 
whereas other authors have taken a priori assumptions which ensure that their 
free boundary stays smooth, we have been able to show existence of weak so- 
lutions for all time, regardless of the possible changes in topology. (Note that 
in |CFj and |FHj they assume that their data is starlike with respect to the 
origin, and in |DL| log concavity of initial data is assumed to guarantee exis- 
tence and smoothness of the solutions.) Indeed there are some very natural 
problems arising in the applications where the topology should change. Con- 
sider for example the problem of what happens with Hele-Shaw flow around 
an obstacle. In this case, there is automatically a change in topology when the 
flow meets itself on the other side. An interesting question is whether or not 
an air bubble will be left behind in the wake of the obstacle. Another obvious 
problem from applications where there will be changes in topology is if there 
are multiple injection slots. In fact, this paper already deals with the second 
situation, since we never assume that D is connected. 

The classical version of the m-approximating problem which we use is given 
as follows: 

1.1 Definition (m-approximating problem). We let u^ m ' denote the so- 
lution of the following partial differential equation, 

«S m) = mA(u {m \x,t) - 1)+, (x,t) ED c x (0, +oo) n {u {m) > 1}, (1.2) 

with boundary data given by 

u M(x,0) = ui(x), x G D c , 
m(u (m) (x,t) - 1)+ = p(x), (x,t) edDx (0, +oo), 

and with the free boundary condition 

(-Vm(ii w -1) + ,(1-H/) )-f = 0, (x,t) E d{vS m) > 1} . (1.4) 

Here v is the outer (n + 1) dimensional normal to the set {u^ > 1}, and this 
free boundary condition will be satisfied when the free boundary is smooth. 
The weak formulation we give at the beginning of the next section will not 
require any regularity assumptions on the free boundary or on the initial data. 
For a fixed m > 0, we call the free boundary problem determined by the 
equations above the m-approximating problem. We will assume 

< ui < 1 is compactly supported, 

< p(x) E C 2 ' a {dD) , (1.5) 

D is a bounded set with 3D e C 2 ' a . 
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1.2 Remark (Simplifying Assumptions). There are two assumptions which 
we make which serve to simplify the exposition but which are absolutely not 
necessary for the derivation of our results: 

1. The assumption that uj has compact support is only technical and serves 
to simplify the construction of certain comparison functions in the fifth 
section. We leave to the reader the verification that it suffices to assume 
only that {uj = 1} is compact for all of the results before the sixth 
section. 

2. The assumption that the Dirichlet data, p(x), on the boundary of the 
slot is a function of x alone can be replaced with the assumption that 
the Dirichlet data p(x, t) is a nondecreasing function of t for each fixed x. 
On the other hand, if p(x, t) is allowed to decrease in time, then Lemma 
(I4.1J1 does not hold in general, and this appears to be essential for the 
methods presented in this paper. 

For convenience, we will define 

M:= \\p\\l°°( 9 d) ■ (1-6) 

Now in terms of the positivity assumption on p(x), we note that p(x) = 
leads to a trivial case: There will be no evolution at all. To see this fact 
extend each by 1 across all of D. In short, the positivity of p(x) is driving 
the evolution. 

In this formulation is energy (or enthalpy), and m(u^ m '(x,t) — 1) + is 
temperature. Equation f 1 1.20 comes from conservation of energy. Basically, 
as m increases, the diffusion happens faster. Competing with this increase in 
diffusion is the fact that the boundary data for on dD is decreasing down 
to 1. As m — > oo we have convergence of our operators to the following picture 
which is typical for Mesa problems: 



4 4 y 




u=0 u=l 



We show that as m — > oo, the converge pointwise to a limit < 
u (°°) < i ; anc l rn(u <jn > — 1) + tend pointwise to a bounded function V(x, t). 
Furthermore, the function V is identically zero in {u^ < 1} and positive and 
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harmonic for fixed t in the component of the set = 1} which contains D. 

(On sets where {u^ = 1} which are isolated from D we will have V = 0, see 
Remark (3.1)1 .) Finally, the pair (u(°°\ V) solves our weak Hele-Shaw problem 
with a mushy region which we give formally in Definition (I2.fi)) . In this way 
we obtain the following results. 

1. We get a natural generalization of the weak formulation of the Hele- 
Shaw problem of DiBenedetto and Friedman (see |DFj ). Our formulation 
allows for a "mushy" region. Moreover, after invoking a result of Bouillet, 
we will be able to say that solutions of our formulation are unique, and 
hence the Mesa limiting procedure we use gives an effective method of 
constructing the solution. 

2. We determine the regularity of the spatial slices of the free boundary by 
invoking the results of Blank, Caffarelli, and Kinderlehrer and Nirenberg 
for the regularity of the free boundary in the obstacle problem (see [Bl , 
P, jKNj). Indeed, if we define 

U(x,t):= fvt\x, s)ds = (l-u I (x))x A(t) (x) (1.7) 
Jo 

where A(t) gives the "puddle" at time t, and 

W(x,t) := f V(x, s) ds , (1.8) 
Jo 

and we formally apply the Baiocchi transformation to the equation 

vt\x, t) = A x V(x, t) , (1.9) 
then for every fixed time to, W satisfies the obstacle problem: 

< W(x, t ) , A x W(x, t ) = U(x, t ) . (1.10) 

In a subsequent paper the authors will address the rectifiability and further 
space-time regularity of the free boundary. 

The paper is arranged as follows: In section 2 we introduce our notion of 
weak solutions to the approximating problem, show some qualitative properties 
of these solutions, derive the existence of the limits (for now in weak-* L°°), 
and give some trivial bounds on these limits. Most of the results in this section 
draw from the maximum principle and from the papers |AK| . |Klj . |K2j . and 
|K3j . Section 3 gives a simple counter-example which motivates the definition 
of the free boundary and of the diffusive region. In section 4 we derive some 
monotonicity properties of our sequences and limits, and as consequences we 
improve our weak-* L°° convergence to pointwise convergence and give an 
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explicit representation of u(°°\ In section 5 we show that the limiting problem 
has a free boundary for all finite time (as opposed to if it "escaped to infinity" in 
zero time or by a fixed time), and we use explicit subsolutions to show that if an 
open ball within the diffusive region has a boundary point on the free boundary, 
then the free boundary will have nonzero velocity at that point. In section 6 
we apply the Baiocchi transformation to Equation (11.9)1 and thereby derive a 
family of obstacle problems that this procedure yields. At that point, under 
suitable assumptions on the initial data, we can invoke the regularity theory 
for the obstacle problem (see jBlj and (Hj) to immediately derive regularity in 
space for the free boundary of our Hele-Shaw problem for almost every time. 
In section 7, we extend the results of section 6 to include every time. In the 
process of extending to every time, we establish the continuity of the measure 
of the diffusive region with respect to time when the initial data is strictly less 
than one. Section 8 is an appendix which gives some of the barrier functions 
and some of their properties that we need in some of the earlier sections. 

I. Blank wishes to thank Kansas State University for their hospitality 
while the bulk of this paper was written. M. Korten is indebted to Juan 
Luis Vazquez for introducing her to the Hele-Shaw problem. M. Korten 
was partially supported by NSF EPSCoR Grant #530517 under agreement 
NSF32169/KAN32170, and I. Blank's visits to Kansas State University were 
also supported by this grant. All three authors wish to thank Barbara Korten 
for her patience and support while this paper was written. 



2 The weak formulations 

For our m-approximating problem, we need an appropriate weak formulation 
which we give here: 

2.1 Definition (Weak solutions of the m-approximating problem). 

The nonnegative function u^ m '(x, t) G Lj oc is a weak solution of the m- 
approximating problem if for any ip G C°° (H n x [0, oo)) , such that ip = 
on dD x [0, oo), and ip(x, t) — > as either t — » oo or \x\ — > oo we have 

ip t (x,t)u^ m) (x,t) dt dx+ / / A x (p(x,t)m[u (m) (x,t) - l] dt dx 
d c Jo Jd c Jo 

= [ f Tr(x,t)p(x) dt dU n ' x x - f <p(x,0)ui(x) dx . 

JdDJo OV J D . 

(2.1) 

We observe that the traces of and m(u^ m ' — 1) + on the boundaries of 
our domain will be well-defined by the work of Korten even if we only assume 
that we are dealing with local solutions and that the initial trace is between 
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zero and one. (See |Klj which adapts the work of Dahlberg and Kenig, |DK| . 
and see Lemma 3.2 of |K2j .) Furthermore, we observe that the free boundary 
condition for the classical formulation (i.e. Equation ( II .4j) ) will be satisfied, 
whenever the functions and sets are sufficiently smooth. 

For strictly local situations, it may be helpful to note that if <p G Qj°(f2) 
where Q CC D c x (0, +00), then we have 



The last integration by parts requires that we invoke the known regularity 
theory and energy estimates for solutions of Equation (11.2)1 . See Lemma 1.2 



Since we want to discuss regularity of functions in Sobolev spaces, we need 
to fix ideas about which representative we will use. Although the "Lebesgue 
point" representative where points are filled in by taking limits of averages 
over balls with radii going to zero is the most common procedure, we will 
use a slightly different approach which exploits the fact that the functions 
u^ m \x, t) are nondecreasing in t for almost every x by Lemma 4.2 of )Klj and 
by the fact that the boundary data p(x, t) is nondecreasing in time for each 
fixed x. 

2.2 Definition (Representative). For t > 0, we set 

u^Hx, t) : = limu (m) (x, s) for a.e. x. 

sit 

For economy of notation we will not bother to relabel any of our w^'s, but 
always understand that we are using as our representative for in all 
pointwise matters. 

We cannot expect to be continuous across the free boundary, and in 
fact, the following theorem summarizes the results of Lemma 2.4 and Lemma 
2.5 of [K3j : 

2.3 Theorem (u^ must jump). The set 




+ 



dt dx 



l] dt dx . 



of UK]. 



E := {x e JR n : 3t s.t. u T (x) < u {m \x, t) < 1} 



has Lebesgue n- dimensional measure equal to zero. 
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On the other hand, m{u^ m > — 1) + will be continuous across the free boundary. 
(See |Klj which verifies the required assumptions of [DB| for continuity.) For 
a fixed m we now have the following picture of viewed from the side: 




Boundary of D 



Free Boundary 



By using the boundedness of the data along with the maximum principle 
that the equation above enjoys, we can conclude that the solutions are bounded 
and therefore by elementary functional analysis we can conclude that the fol- 
lowing limit 



u 



(oo) 



lim u {m \x,t), 



(2.2) 



exists weak-* L°° along a subsequence of m in the entire domain. Again by 
using the maximum principle it is immediate that 



< u {oo) < 1 . 

In fact, by using the maximum principle again, we can assert that 



< u <«0 < 1 + - 
m 



or 



< m(u {m) - I) , < M 



(2.3) 



(2.4) 



and we stress that M is independent of m and t. Because of this fact, we can 
take a further subsequence to ensure that m(u^ m ' — 1) + has a limit V in the 
weak-* L°° topology of the entire domain. 

2.4 Remark. By combining Equation (12. 4|) with Theorem (|2.3J1 we can con- 
clude that the essential range ofu^ m '(x, •) is a subset of {itj(x)}U[l, 1+M/m] . 

Next, for (p G Hq(Q) we have: 



/ / V, [m{u {m) - 1)+] • dxdt= / 

JO J D c JO J D c 



JD 



u {m) Vt dx dt 
w (oo Vt dx dt . 



Since u(°°) G L°°(fi) C L 2 (Q), we have 

u ( r ] g H-\n) 



(2.5) 
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We summarize with the following lemma whose proof is now obtained trivially 
by using the weak-* L°° compactness we have and taking the limits as m — > oo 
in the weak formulation of the m- approximating problem. 

2.5 Lemma (Limits solve the limiting problem). Under our assumptions 
as above, we have 

< V < M , (2.6) 

and the pair (u^, V) satisfies 

p poo 

(pt(x,t)u ( - 00 \x,t) dt dx + / / A x ip(x,t)V(x,t) dt dx 
d c Jo Jd c Jo 

(2-7) 

= f [ *rf-(x,t)p(x) dt d^^x — [ (f(x,0)uj(x) dx . 

JdD JO UV J D c 

for any if e C°° (H n x [0, oo)) , such that if = on dD x [0, oo), and ip(x, t) — > 
as either t — > oo or \x\ — > oo. 

2.6 Definition (Hele-Shaw with mushy region). Any pair (u^°°\ V) 
which satisfies Equation ( 12. 7j) for if as in the lemma above will be called a 
weak solution of the Hele-Shaw problem with boundary data p(x) and initial 
data ui{x). 

2.7 Remark (Inclusion of earlier models). The weak formulation we have 
above generalizes the Hele-Shaw formulation of DiBenedetto and Friedman to 
include the situation where there is a mushy region (see |DF| ). In addition, our 
a priori regularity assumptions on the solutions are weaker. (We only assume 

2.8 Theorem (Uniqueness among all solutions). The solution of the 
limiting problem of Lemma H2.d}) is unique, so all solutions of the Hele-Shaw 
problem as given in Definition \2. 6}) are recoverable via the Mesa limit process 
we have introduced. 

Proof. This theorem is an immediate application of the uniqueness theorem 

in section 3 of |Boj . 

Q.E.D. 



3 Counter-example to regularity in time 

Based on the assumptions we have so far, the function V(x, t) will not be 
continuous in time in general. To show this, we will assume for the sake of 
this example that ui is a continuous function. We will also assume in this 
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section that the set D is connected. The trouble appears in certain cases 
where the set 

W := D U {x E H n : uj(ar) = 1} (3.1) 

is disconnected. If K is a component of the set {uj = 1} and K is a positive 
distance away from D, then the and therefore the should not evolve 
on this set of x until the free boundary comes into contact with it. Essentially, 
the "patch" K will not "see" the input of the slot until the component of the 
set {u^ > 1} which surrounds D connects to it. 

3.1 Remark (V is not always continuous). To produce a situation where 
V must be discontinuous in time, simply consider the radially symmetric sit- 
uation where D = B\ and ur(r) is taken such that it is identically one on the 
set 3 < r < 5, but smaller outside of it. When what we want to call the free 
boundary reaches r = 3, then it will instantaneously jump to r = 5, and this 
leads to an immediate jump in the height of V. 



The free boundary at different times. 




Assume that the initial 
data is exactly one in 
this set, and less than 
one outside of it. 



Before this contact, V is identically zero on 
all of K. 



With these pictures in mind we make the following definitions: 
3.2 Definition (Defining the free boundary). We set 

t rn (x) : = inf{t : m(V m )(:r,t) - 1)+ > 0} 
too (a;) := inf{* : V (x,t) > 0} 

and 

T m (x) := inf{f : w ( "%,t) > 1} 
^(a;) := inf{t : u^(x,t) = 1} . 



(3.2) 



(3.3) 



We can also now define the diffusive region at time t for the m-approximating 
problem and the limiting Hele-Shaw problem to be respectively, 

A^ m \t) := {x G D c : m(u {m) (x, t) - 1)+ > 0} , . 

A(t) := {x E D c : V(x, t) > 0} . [ ] 
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Now we simply define the free boundary at time t to be FB^ (t) := dA^ m ' (t) \ 
dD, or FB(t) := dA(t) \ 3D, according to which problem we are considering. 
(Here u dK" denotes the topological boundary of K.) 

Observe that because of the radially symmetric example given above where 
V is discontinuous in time, it is also clear that the free boundary cannot be 
expected to vary continuously in time. 

3.3 Remark ("Diffusive" versus "Wet"). The reader may observe that 
"Wet" may be a more appropriate name for what we are calling the diffusive 
region in the Hele-Shaw problem. On the other hand, "wet" might more 
properly describe any region where t) > 0, so we will stick to what is 

already appropriate in the approximating problems. 

By considering the figure above, it is clear that in general t m (x) does not 
have to equal r m (x), and similarly with t^x) and Too(x). In the figure, they 
would differ on the set K. For x G K we have r m (x) = t^x) = 0, while 
t m (x) would be the positive time when the moving part of the boundary of 
| M (m) > crosses K, and t^x) would be the positive time when the moving 
part of the boundary of {u^ = 1} crosses K. 

4 Monotonicity, inclusions, and consequences 

We start with a lemma which is a simple consequence of Lemma 4.2 in [KlJ. 

4.1 Lemma (Monotonicity of u^). For each m > and any x e D c we 

have vS m \x,t) is an increasing function of time. 

Proof. If (x, t) < 1, then this follows immediately from Lemma 4.2 of |Klj . 
As a consequence of this fact, we can already say that the A^ m \t) are nested 
nondecreasing sets. Now we consider the set S := {(x, t) : U ( m \x,t) > 1}. By 
observing that satisfies the heat equation within S and invoking standard 
parabolic regularity theory, we can conclude that t) will assume the 

boundary values 1 + on dD continuously, except possibly at the initial 
corner. By Corollary 1.3 of |AK| the function («("*)- 1)+ is continuous within 
D c , so we can extend to be 1 on the rest of the parabolic boundary of S 
in a continuous fashion, except again at t = where this set meets dD. Now 
we consider the function 

U{x,t) :=u {m \x, t + e) -u {m \x, t) . (4.1) 

Within S, this function satisfies the heat equation. On dS fl dD, we have 
U(x,t) = 0. (If we allow p(x) to also depend on time, but we require it to 
be nondecreasing in time, then we have U(x, t) > here which is also fine. 
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It is not at all clear, however, how to deal with the most general case where 
p(x,t) is allowed to decrease in time.) On the rest of the parabolic boundary 
of S we have U(x,t) = u^ m '(x, t + e) — 1 > 0. In the corner where U(x,t) is 
discontinuous, it stays bounded between zero and 1 + and hence by the 
weak maximum principle, U(x,t) is nonnegative in all of D. Hence 

u^(x, t + e)-u^ m \x, t) 

~ (4-2) 

is nonnegative, and the result follows by taking e | 0. 
Q.E.D. 




4.2 Lemma (Monotonicity of u^). is monotone increasing in t for 
C n a.e. x. 

Proof. By the last lemma the are monotone increasing in t. Now to 
come to a contradiction, suppose h > is a number and (fi(x,t) is a bounded 
nonnegative function which is compactly supported in our domain and which 
satisfies 

'**,«)„<-><*,«)* *>//„<*, «-%<-><,, t)dxit . 

Such a construction is possible if the conclusions of our theorem are not satis- 
fied. Then ip(x,t) := ip(x, t) — (p(x, t — h) is an admissible test function which 
satisfies 

//«,,*)«<">(*,*)**><., and //«*,«)„<->(*, t)(fa *<0 

which contradicts the weak-* L°° convergence of the to u^. 
Q.E.D. 

Now without data on the slot which "competes" with the increasing diffu- 
sivity (or if the functions solve the equation on all of ]R n ) , it is immediate by 
rescaling that the diffusive regions (the A^ m \t)) must be nondecreasing with 
m. In our case, however, we need to produce an appropriate barrier. 

4.3 Theorem (Temperature increases with m.). Ifm<k, then 

m[u (m) - 1]+ < k[u (k) - 1]+ . (4.3) 

As a trivial consequence we can say that the A (m \t) are nested: m < k implies 
A^ m \t) c . 
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Proof. We make the following definition: 



( u( m \x, t) if u (m \x, t) < 1 

vW(x, t) := I (4.4) 
{ 1 + f («("%, t) - 1) if i) > 1 

and claim that v^ k \x, t) is a subsolution of the k- approximating problem. 
Indeed, this follows very quickly from the following two easily verifiable state- 
ments: 

1. The set where v ^ > 1 is identical to the set where > 1. 

2. k{v^ - 1)+ = m(u {m ^ - 1)+. 

(Only the definition of v^ k ' is needed to verify these statements!) By the first 
observation the free boundaries of v^ k ' and are identical in time and space, 
so the speed of these boundaries at every point is also identical. On the other 
hand, the second observation quickly leads to the conclusion that the function 
v ^ satisfies the free boundary condition and the boundary value condition on 
the slot for the k- approximating problem exactly. Finally, we simply compute: 

v?\x, t)< M J m) (x, t) 
= mA{u im) {x, 
= kA{v {k \x, t) - 1)+ 

which shows that v ^ is a local subsolution to the k-approximating problem, 
and therefore v <yk \x, t) < u^(x, t). (In the first inequality above we have 
used Lemma (j4.1j) .) We now have the simple consequence 

k[u^ - 1]+ > k[vW - 1]+ = m[u^ - 1] + , (4.5) 

which is the monotonicity we require. 
Q.E.D. 

4.4 Corollary (Pointwise convergence of the temperature). The se- 
quence of functions {m[u^ — 1]+} converges pointwise almost everywhere to 
V. In particular, the limiting function V is unique. (No subsequence is ever 
needed.) 

Proof. By the last theorem combined with the estimates in place already, at 

each point we have a bounded increasing sequence. 

Q.E.D. 
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4.5 Corollary (Representation of u^). There is an increasing set-valued 
function oft which we call Q(t) such that t) admits the representation 
for almost every (x, t): 

u {oo) (x,t) = X Q{t) (x) + ui(x)x Q(t)C (x) . (4.6) 

Furthermore, t) is the pointwise limit of the functions u {m \x, t) al- 

most everywhere, and Q(t) can be chosen to be equal to the set {x G R n : 
t~oo(x) < t] (Q{t) is increasing, means in the sense of set inclusion.) 

Proof. By Lemma (14.2)1 . we know that u^°°\x, t) is an increasing function 
of t and so by using Lemma (J4.2j) along with Equation (12. Hj) we can conclude 
that ui{x) < u^°°'(x, t) < 1. So, to prove Equation ( 14 .6 j) it suffices to show 
that ■) does not attain values strictly between uj(x) and 1 for a.e. x. 

Indeed, by combining the previous theorem with Remark 1)2 .4j) we see that the 
limit 

lim u^(x , t ) (4.7) 

m— >oo 

exists for almost every (xo,t ) G H n+1 and is either equal to 1 or to Ui(xo)- 
At that point it is a simple exercise to show that this pointwise limit coincides 
with the weak-* L°° limit almost everywhere. 

In terms of the "choice" of Q(t), it is clear that the only x where there can 
be a choice is on the (possibly empty) set {x G IR n : Ui(x) — 1}. Examina- 
tion of the definition of Q(t) combined with the monotonicity of u^°°\x, t) in 
t makes it clear that := {i 6 R n : t^x) < t} will suffice. (Note that an 
equally good choice for Q(t) would be to take Q(t) := {x G R n : t^x) < t}.) 
Q.E.D. 

4.6 Remark (Uniqueness of limits of subsequences). Although we had 
originally needed subsequences to be sure that our would converge in 
weak-* L°°, the last result makes it clear that will converge both pointwise 
and weak-* L°° to without the need to extract a subsequence. 

4.7 Corollary (V is harmonic within the diffusive region). The spatial 
Laplacian of V is zero for x in the interior of A(t). As a consequence, if 
3D G C k ' a where k > 2, then V attaches to the slot D in a C k,a fashion. 

Although the function V is discontinuous in general, we do get regularity in 
space. 

4.8 Theorem (Spatial continuity of V.). V( ■ , t) is continuous for almost 
every time, t. 
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Proof. We will prove that V( ■ , t) is continuous by proving that it is both 
upper and lower semicontinuous. Lower semicontinuity follows immediately 
from the fact that the functions m[u^ m \ ■ , t) — 1] + are continuous and are 
increasing as functions of m. To get upper semicontinuity we have a little bit 
more work. 

By Lemma 1)4. 2|) . dv = u[°°^ is a nonnegative measure. Our Radon measure 
can be "sliced" into Radon measures of one dimension less for a.e. t. (See [M] 
p. 139 - 142, and Equation (10.3) on p. 140 in particular.) We will call these 
slices du t . Since for a.e. t we now have A x V(x, t) = du t , we know V(-,t) is 
subharmonic, and therefore upper semicontinuous for those values of t. 
Q.E.D. 



5 Nontriviality and nondegeneracy 

Due to the competition between the decreasing data on 3D and the increasing 
diffusivity of our m-approximating problems, we need now to rule out two 
trivial cases. The possibilities that we need to exclude are: 

1. The possibility that the free boundary moves to infinity as soon as t > 0. 

2. The possibility that the free boundary never moves. 

Conceptually, our elimination of these cases is trivial. We simply produce 
a family of supersolutions (and then subsolutions) to our m-approximating 
problems (Equation (jl.2j) ) whose free boundaries move in a suitable manner, 
independent of m, and then we make use of the comparison principles that 
our equations enjoy (see |Klj ). 

5.1 Theorem (Boundedness of the free boundary). If ui is compactly 
supported, then so is u^°°\x,to) for any to > 0. In fact, no matter what the 
initial data, there is always a supersolution whose free boundary has constant 
speed. 

Note that we need the compact support of the initial data for this result to 
hold. 

Proof. By using the boundedness of D, and by translating and rescaling 
it will suffice to produce a supersolution to the problem where D c = B\, 
ui = x { b 2 \ Bi }' an d p( x ) = k- We claim that the function 

- ,m) (^) : K^(Sr +1 )^2 + «] w (51) 

will suffice if I is sufficiently large. Note that this function is linear in the 
radial variable, and that the free boundary moves with speed i for all time. 
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The following computations are elementary, and they prove our claim that 
v ( m > is a local supersolution: 



m 



A(v {m) - 1). 



k(l-n) k£(r - 1) = (m) 

r(l + £t)~ -m(l + £t) 2 1 



(5.2) 



If we are on the free boundary, so that r = 2 + it, then 







- [m{v™ - 1) + ] 



[l + £t) 



< £ = speed of the free boundary. (5.3) 



Q.E.D. 



5.2 Theorem (Nontrivial motion of the free boundary). The boundary 
of the diffusive region does not remain stationary, and for any R > 0, there 
exists a time t such that for all t > t , we have Br C A(to). 

5.3 Remark. If we only want to show that the free boundary does not remain 
stationary, then we can use the same subsolutions as in the proof of Theorem 
(|4.3|) . On the other hand, the subsolutions we use here are also needed in the 
proof of the next theorem. 

Proof. Using the fact that D contains an open set, and by translating and 
rescaling it will suffice to produce a subsolution to the problem where (for any 
a > 0,) D c = Bi +a , Ui = Y rn . „ and p(x) = k, and such that the free 
boundary of this subsolution moves a fixed distance (independent of a) in a 
finite time T(a). 

We will create our subsolution in a couple of steps. First, let w^ m \r,t) = 
wf 71 ^ (r, t) + u>2 ( r ? t) where (r, t) solves for each fixed t 

A x w4 m) (r, t) = in B 2+a+u \ B 1+a 

wt\r,t) =1 + - ondB 1+a (5.4) 
m 

u;S m) (r,t) =1 ondB 2+a+ et 
and (r, t) solves (again for each fixed t) 

in B<2+ a +tt \ B\ +a 

(5.5) 

on d{B 2+a+ £t \ B 1+a } . 



A x w { 2 m \r,t) =- 
m 

w [ 2 m \r,t) =0 
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Since the free boundary will be given by r = 2 + a + it, and we only need it 
to move a fixed distance, we will also assume that it < 1. The functions w[ m ^ 
and u>2 can t> e gi yen explicitly, and their relevant properties are given in the 
appendix. 

By invoking Corollary ()8.3|) from the appendix, we have 

9 

-oo < -C x {ri)k < —mw ( { m \2 + a + it, t) < -C 2 {n)k < (5.6) 

where we stress that the constants are independent of a and it. (See the 
appendix for a few more details, and remember that < It < 1.) By the same 
corollary we can conclude that 

r) 

< C 3 (n)e < —mwt ] {2 + a + it, t) < CAn)e < oo (5.7) 
or 

again with constants independent of a and it. Now we take e > sufficiently 
small to ensure that 

-oo < -dk < ^- \mw {m) {2 + a + it, t)] < -C 2 k < . (5.8) 
or L J 

Now by taking i < Cikjl we can be sure that our function is a subsolution 
along the free boundary. We note that Am(w' m '(r, t) — 1) + = e > in the 
region where w^ m \r, t) > 1. Since 

lim wl m \r, = 0, (5.9) 

m— >oo 

Once m is sufficiently large, we automatically have 

Am(w {m) (r, t) - 1)+ > w t (m) (r, t) . (5.10) 

Since the free boundary of our subsolution moves with speed i > 0, we are 

done. 

Q.E.D. 

The following theorem shows the instantaneous detachment of the free 
boundary from the slot, dD, even if uj(x) < A < 1 in all of D c . 

5.4 Theorem (Instantaneous formation of the diffusive region). // 

to > 0, The set A(t ) contains an open neighborhood of dD. 

Proof. Because we have assumed that dD G C 2 ' a , every point on dD can 
be touched from within with a tangent ball, and then we can use the same 
subsolutions of the previous theorem to force instantaneous movement of the 
free boundary. 
Q.E.D. 
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6 Spatial regularity results 

In this section we will derive spatial regularity for both the function V(x,t) 
and for the free boundary. We apply the Baiocchi transformation to V(x,t) 
and define: 

W{x,t) := [ V{x, s) ds . (6.1) 
Jo 

Observe that by Lemma ( 14. 2 j) (which shows that the diffusive regions are 
increasing in time) and by the positivity of V(x,t) in the diffusive region, it 
is clear that the set {W > 0} is identical to the set {V > 0}. Now to find the 
regularity of d{x G !R n : W(x,T) > 0} we will show that W(-,t) belongs to 
H} oc (D c ) for almost every t, and then that W(x, t) is a weak solution of the 
following obstacle problem in almost every time slice t = T 

0<W{x,T), A x W{x,T) = x, m 

x,t)>o} ( x ) (1 — u i( x )) ■ (6-2) 

After that we will be able to invoke regularity results for the obstacle problem 
due to Caffarelli, Kinderlehrer, Nirenberg, and Blank. 

For simplicity, we let a m (s) := m(s — 1) + . We start by stating some simple 
trace results. Basically, we need to adapt Equation (12.1(1 to some situations 
with slightly different test functions. 

6.1 Lemma (First Trace result). If ip{x) G C°° is supported in the interior 
of D c , then the following formula holds for a.e. T: 

[ [ [A x ip(x)) a m (u^ m \x,t)) dt dx 

J D c JO 

= [ ^(x)[w (m) (x,T) -u^x)] dx . 
Jd c 

Proof. We make the following definition 

( 4){x) t < T 
<p{x,t):=\ (6.4) 

[ e(x,t) t > t 

where Q(x,t) is chosen to ensure that (p(x,t) is a permissible test function for 
our m-approximating problem. In particular, we need ip(x,t) G C°° and we 
need it to converge to zero as t — ► oo. Neither requirement poses any difficulty. 
By using the trace result of jAKj (see Theorem 1.1 of |AK| ) our functions 
t) solve our m-approximating problem starting at time T with initial 
data u^ m '(x,T) for almost every T, and so (for those T) we can use Q(x,t) as 



(6.3) 
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the test function in Equation (|2.1|) to obtain: 

Q t (x,t)u {m) (x,t) dt dx 



D c JT 



D c JT 



oc 



+ / / A x e(x,t)a m (u {m \x,t)) dtdx (6.5) 



i){x)u {m \x,T) dx . 

(Note that Q(x,T) = i/j(x) since we have required that <p(x,t) be smooth.) 
Now by subtracting this equation from what we have when we plug in the 
function ip(x,t) defined in Equation l|6.4j) into Equation (|2.1j) we get Equation 
()6.3|) immediately. 
Q.E.D. 

The proof of the following result is almost identical to the proof of the 
trace result above, so we omit it. 

6.2 Lemma (Second Trace result). For a.e. to,h such that < t < 
t\ < oo, and for ip 6 C°° (!R n x [t ,£i]) , which satisfies (p = on an open set 
containing D x [io^i], we have 

h r rti 

ip t (x,t)u {m \x } t) dt dx + / / A^(:r,t)a m (w (m) (:r,t)) dt dx 



D c Jt JD C J to 

ti s 



[ V (x,s)u {m) (x,s)] ) dx. 

s=t J 

(6.6) 

Now we state the standard energy estimate for our situation. 

6.3 Lemma (Energy estimates for the uS™^). Let < r < R and < t < 

t\. Then there is a constant of the form 

C(n) , . 

° = jR—rf (6 - 7) 

such that the following energy estimate holds: 

I 1 I \Va m {u {m) )\ 2 dx dt < C [ 1 [ a m {u (m) ) 2 dx dt . (6.8) 

J to JB r {x ) J t JB R (x ) 

6.4 Remark (Independence of Time). Notice that the constant C is in- 
dependent of time and notice that the time intervals in the integrals in each 
side of the inequality are identical. 
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Proof. We choose rj(x) £ C£°(B R ) such that 

1. r) = 1 on -B r , 

2. < 77 < 1, and 

3. |Vr/| < 4( J R-r)~ 1 . 

Now let ip(x,t) := a m (u^ m ')r](x) 2 , and apply the last lemma and Green's iden- 



tity to obtain: 

<Pt(x, t)-u ( - m - ) (x, t) dt dx — 

B R Jto J B R 



[(f{x,s)u im) {x,s)} 



s=ti 



s=to 



dx 



(6.9) 



Br Jt 



S/ X ip(x, t)V :r a m (-u (m) ) dt dx . 



The fact that we can apply Green's identity above can be justified using Lemma 
1.2 of |AKj . By Lemma lj4.1J) we know that for a.e. x £ D c , fi x (t) := Ut(x, t) is 
a non-negative Radon measure, whence it suffices to have <f>(t) £ C (0, +00) for 
the distributional pairing (fJL x (t), <j>(t)) to be defined, and non-negative if > 0. 
Using Equation j2.1J) it is easy to see that the function g(x) := (fJ- x (t), X[ t H ] (0) 
is locally integrable in x. In other words, we can integrate the left hand side 
of Equation ()6.9|) by parts in time to give us the following inequality: 



Br Jt 



ti 



V, (a m (u (m V(a:)) V x a m (u (m) ) dt dx 



(a m {u^)r, 2 {x))d^ x {t) 



to 



dx 



< 0. 

This inequality implies 



Br Jto 



n 



2 : |V \a m (u {m) )\ 2 dtdx < 



1 I (v IV^u^)!) {a m (u {m) ) \V xV \) dt dx 

'Br Jto 



and so Cauchy-Schwarz gives 
'■ti 



to 



< 4 



V( x ) 2 |V x a m (w (m) )| dx 



UB R 
ti 

to 



dt 



(6.10) 



b. 



\V xV (x)\ 2 a m (u {m) ) 2 dx 
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dt . 



Q.E.D. 



The aforementioned independence of time of the constant C leads imme- 
diately to the following theorem: 

6.5 Theorem (Time-slice energy estimates). For almost every time t we 
have 

[ \Va m {u {m) {x,t))\ 2 dx < C [ a m {u {m) {x,t)f dx , (6.11) 

JB r (xo) Jb R {x ) 

and also 

f \VV(x,t)\ 2 dx <C f V(x,t) 2 dx, (6.12) 

J B r (x ) J B R (x ) 

where once again, the constant C has the form given in Equation \6. 7\) . 

Proof. Equation ( I6.11j) is already known from the previous lemma, and 
Equation (J6.12J) follows from Equation . 1 lj) by using the L 2 convergence of 
®m(u( m \x, t)) to V(x, t) and by using the lower semicontinuity of the Dirichlet 
integral. 
Q.E.D. 

6.6 Corollary (Energy estimates for V(-,t)). For almost every t, V(-,t) G 

6.7 Theorem (Energy estimates for W(-,t)). For all t > 0, W(-,t) G 
HL(D C ). 

Proof. Let K be a compact subset of D c . Since W(-,t) < Mt, we have 
W(-,t) G L°°(K) C L 2 (K). It remains to show that f^(-,t) G L 2 (K). 
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Let ip G Cq(K) and estimate. 



dW 

dxi ' 



K 




W(x,t)^-(x) dx 



V(x, s)^-(x) ds dx 



K JO 



dxi 




k Jo 



— lim 



lim 



lim a m {u (m '{x, s)) — — (x) ds dx 

TTW + OO ' OXi 

am{u < - m \x, s)) tt^-(x) ds dx 




A JO 



dxi 




1 d 



'A JO ^Xi 



{a m (vA m \x } s))) (p(x) ds dx 



lim / <p(x) 

m— >+oo 



A' 



J A {a m (u (m) (x, s))) ds 



dx 



where we have used Lebesgue's Dominated Convergence Theorem repeatedly 
in the computation above. Now, by using Minkowski's integral inequality and 
the Cauchy-Schwarz Inequality we have 



<p(x) 



A 



d_ 

o 9xj 



{a m (u (m \x, s))) ds 



dx 







1/2 




< 


/ (p(x) 2 dx 




Ik ^kI 




JK 







jn^ {a m {u^\x } s))) ds^j dx 



1/2 



< 



< 



m\L 2 (K) 



L\K) 



A 



-^L(a m (u {m \x, s))) ) dx 



1/2 



JK 



{oL m {u^ m \x } s))) ) dx ds 



ds 




1/2 






/ ds 




Jo 



1/2 



< t 1/2 \\Lp\\ L 2 {K) \\\7 x a m (u {m) (x, s))||i2(Kx(o,t)) 

< Ct 1/2 \\^\\ L 2^ 

where the last constant is independent of m by Lemma ()6.3j) . 
Q.E.D. 



6.8 Remark. In fact, once we show that W(-,t) is a solution of the obstacle 
problem, we will be able to infer from elliptic regularity theory, that W(-,t) is 
C l,a in space for all a < 1. 
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In order to derive Equation (16. 2|) we will need to commute the Laplacian 
with the integral in time. To accomplish this commutation we turn back to 
the approximating problem where this switch is simpler. 

By using weak-* L°° convergence of the temperature functions we have 



J D c 

On the other hand, by using the dominated convergence theorem we conclude 
that for almost every T 



Since W > 0, we can combine the last two computations with the previous 
lemma to conclude that W(-, t) solves the obstacle problem for a.e. t > 
0. Equation ( 16. 2|) and this fact allows us to use the technology of [B^ and 
[H] to infer regularity of the free boundary here as long as we satisfy the 
nondegeneracy condition: 



(ui satisfying Equation ^6.13j) will be referred to as nondegenerate initial data.) 
In this case we will have the next theorem which we state after one simple 
definition. 

6.9 Definition (Minimum diameter). The minimum diameter of a set 
S C JR n (denoted u m.d.(S) v ) is the infimum of the distances between parallel 
hyperplanes enclosing S. 

6.10 Theorem (Regularity of the free boundary in space). Assume uj is 
continuous, nondegenerate initial data. Then there is a modulus of continuity 
a which depends on A, n, and the modulus of continuity of uj such that at 
almost any time t, (indeed every time t where Equation \6 . fy) is valid) and for 
any free boundary point (x ,to) contained in the interior of D c we have either 
(with B r (x) denoting the (spatial) ball centered at x with radius r) 





ui(x) < A < 1 . 



(6.13) 




for all r < 1 , 



(6.14) 



23 



(so that the nondiffusive region is "cusp-like") or 

r^o \B r \ 2 v ; 

in which case we will say that (xq, to) is a regular point of the free boundary. 

Proof. Simply use Equation ()6.2j) together with the results from [Blj . 
Q.E.D. 

6.11 Theorem (Better regularity). Assume the hypotheses of the previous 
theorem and assume that (xq, to) is a regular point of the free boundary. Then 
the free boundary intersected with {t = to} will be Reifenberg vanishing nearxo- 
(For a definition of Reifenberg vanishing, see {Bl}.) If uj is Dini continuous, 
then the free boundary will be C 1 (in space) near xq, and if ui G C k,a then the 
free boundary will be C k+1 ' a (in space) near xo- 

Proof. Again, just combine Equation (j6.2j) with the results from |Blj . 
Q.E.D. 

6.12 Remark (Clarification). When we say "near" x in the theorem above 
we mean near in space only. In other words, we are specifically talking about 
the free boundary restricted to the time slice t = t . 

At this point we need to call attention to our assumptions that uj is con- 
tinuous and nondegenerate. Indeed it is well known that there are examples 
of "persistent corners" in Hele-Shaw problems. (See |KLVj for example.) Our 
theorem does not contradict this fact. Our spatial regularity theorem says 
nothing if uj is discontinuous, or if the set {ui = 1} C\D C is nonempty. On the 
other hand, our assumptions do not rule out "focusing" or changes in topology 
of the diffusive region, so it is certainly nontrivial that corners do not arise 
in this setting. Consider where ui is very close to zero except for an e 

neighborhood of an annulus which contains the slot. In this annulus assume 
that Ui is extremely close to one. (The epsilon neighborhood is needed to make 
ui continuous.) In this case, the diffusive region will make its way around the 
annulus in each direction very quickly, but expand slowly into the region inside 
and outside of the annulus. Thus, it will meet itself on the other side of the 
annulus long before it fills in the interior. 
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7 Continuity and dealing with all times 

We wish now to improve the results from the previous section by extending 
them from almost every time to every time. Indeed, this prevents corners from 
arising as "transitions" from the situations which are permissible for sets of 
time with positive measure. (Consider for example the sets S(t) := {(x,y) : 
xy < t — 1 for the interval t G [0, 2]. By the results of the last section, it would 
still be possible for S(t) to be the diffusive region for a Hele-Shaw flow, since 
the corner only occurs when t belongs to the zero measure set {1}.) Of course in 
this entire section we make the standing assumption that ui is nondegenerate 
initial data (see Equation Ijfi.lHjl ). We start with a simple lemma summarizing 
some of the regularity we have for W. 

7.1 Lemma (Regularity for W). W(x,t) is continuous in space and con- 
tinuous and convex and nondecreasing in time. For almost every time, T, 
W(x,T) satisfies Equation \6. fy) with boundary data 

W(x,T)=p(x)T xedD. (7.1) 

Proof. By using the maximum principle together with the fact that the dif- 
fusive region increases with time, we see that V(x, t) must be an increasing 
function of time. Now convexity of W follows from this fact and from the 
definition of W. The rest of the lemma follows immediately from the definition 
of W and the spatial continuity of V(x, t) along with its boundary data on 
dD. 
Q.E.D. 
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7.2 Theorem (Measure of the diffusive region). There exists a constant 
C = C(n, a, 3D) such that ifO < t — s < 1, then 



\A(t)\A(s)\<C\\p\\ c ^ {dD) lt_^ 



(7.2) 



(Recall that A(t) is the diffusive region at time t, and for S C IR", we let \S\ 
denote the Lebesgue n- dimensional measure of S.) 



Proof. We adapt the proof of Theorem 4.1 of [Bl to the current setting. Fix 
t and s. Because the diffusive regions are nested, and because Equation (J6.2j) 
holds for almost every time, without loss of generality we can assume that it 
holds for both t and s. Let L := A{t) \ A(s), let := W(x,t) - W(x,s) 

and observe that A$f(x) = x L (l — ui) > 0, and 



= 16 FB(t), 

qr(x) = (t - s)p(x) x e 3D . 



(7.3) 



From this fact and by the weak maximum principle, it follows that 

< W(x,t) - W{x,s) < (t-s)\\p\\ L ~ (dD) (7.4) 

for all x G A(t). (Nonnegativity is actually a consequence of the previous 
lemma. ) 

Observe that W(x,t) — W(x,s) is harmonic within A(s) so that 



0= / A(W(x,t) - W(x,s))dx 

lA(s) 

d_ 

i;>D dv 

d_ 

i; id < )ly 



(W{x,t) - W{x,s))am n - 1 
(W(x,t) - Wix^^dH 71 - 1 



FB(s) 



FB(s) 



d_ 

dv 
d_ 
dv 



(W(x,t) - W(x, s))dH 



n-1 



W(x,t)dH 



n-1 



Now by using boundary regularity for harmonic functions combined with Equa- 
tion (|7.4|) we can conclude 







(W(x,t) - W(x,s))dH 



n-1 



< C(n,a,dD)(t-s)\\p\\ 2„ 



idD dv 

By combining this fact with the last computation, we conclude that 
d 



(8D) 



(7.5) 



/ 



dv 



W(x,t)dH 



n-1 



< C(n,a,dD)(t - s)\\p\\ C 2 H 8D) ■ (7.6) 
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On the other hand we have 

;i-A)|L| < / (1- Ul )dx 

AW(x,t) dx 

^-W(x,t)dH n - 1 
dL ov 



[ -^-w(x,t)dn r ' 

Jfb(s) du 



in— l 

7^— vv yx, i jar 

1(8) 

which we can combine with Equation ()7.6|) to give us what we need. 
Q.E.D. 



7.3 Corollary (Continuity in L p ). Under the assumptions made at the 
beginning of this section, the map from t to the function 

X{w(x,t)>0} 0- ~ Ul ) 

is a continuous function from IR into L P (D C ) for 1 < p < oo. 

7.4 Corollary (Spatial regularity for every time). All of the results of 
Theorems l\6'.l(J\) and ifj.ll\) hold for every time. 

Proof. It suffices to show that Equation ( 16. 2j) holds for all time. Fix t > 
and let t n — > t with t n chosen so that Equation f l6.2j) holds at each t n . Now 
take a ball, Br which is large enough to contain A(t) in its interior, and let 
Q := Br \ D. For each n we let w n (x) solve the boundary value problem 

Aw n (x) = x A(tB) (l- wj) in Q 

w n (x) = t n p(x) on dD (7.7) 

w n (x) = on 8Br . 

By standard uniqueness results, w n (x) = W(x,t n ), as they satisfy the same 
boundary value problem. By standard elliptic regularity theory, since the 
boundary data on dD will converge to tp(x) and by the last corollary the 
right hand side of the equation will converge in L P (Q) to X A(t -)(l — u i)i we can 
conclude that w n will converge to a function w which satisfies 

Aw(x) = X A(t - } (l - «/) in n 

w(x) = tp{x) on dD (7-8) 

w(x) = on 8Br . 

On the other hand, w n (x) = W(x,t n ) converges to W(x,t) by the continuity 

of W in time. 

Q.E.D. 
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8 Appendix 

Here we collect some facts we need to construct our subsolutions in the proof 
of Theorem (|5.2j) . We let u(r; a, (3) denote the solution to: 

A x u(r;a,f3) =0 in B 1+a+/3 \ B a 
u(r;a,P) =1 on dB a 
u(r;a,j3) =0 on dBi +a+/3 
and we let v(r; a, (3) denote the solution to: 

Av(r;a,P) = 2n in B 1+a+p \ B a 

v(r;a,P) =0 on d{B 1+a+p \ B a } . 
We will always assume that a > 1 and < (3 < 1. 

8.1 Lemma (Explicit Forms of Our Comparison Functions). If n > 2, 

then u has the explicit form 



(8.2) 



2-n 



u(r; a, (3) 



and v has the explicit form 



v(r; a, (3) = (r 2 — a 2 ) + 



r 2 ~ n - (1 + a + (3) 

a 2-n _ (J + ~ + pj2-n 



{a 2 - (I + a + (3) 2 )(r 2 - n - a 2 - n ) 
(1 + a + (3) 2 - n - a 2 ~ n 



(8.3) 



U) 



If n = 2, then u has the explicit form 



u(r; a, (3) 



log 



1 + CK + /3 



log 



1 + a + (3 



(8.5) 



and u has the explicit form 



v(r;a,/3) = (r 2 - a 2 ) + 



log (£) (a 2 -(l + a + {3f 



log 



1 + a + (3 



a 



(8.6) 



8.2 Lemma (Derivatives on the Outer Boundaries). Ifn>2, then 

(2-n)(l + a + /3) 1 -™ 
w r (l + a + (3- a, /?) = \ — ■ ^[ 2 _ n < , 



a 2-n _ (1 _|_ a + ft} 



1.7) 
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and 

v r (l+a+(3; a, /3) = 2(l+a+(3) + (n-2)(l+a+(3) ^±^±^ ~ ^- > . 

(8.8) 

If n = 2, then 

u r (l + a + P;a,P) = l —, ^<0, (8.9) 



(1 + a + P) log 

and 



a 



1 + a + (3 



u r (l + a + /3;a,/3) = 2(l + a + /3) + ~ l ,, P) > . (8.10) 

(l + a + P) log ' 



a 

Also, for any n > 2, we have 

lim u r (l + a + /3;a,P) = ^<0, (8.11) 

a^oo 1 + p 

and 

lim u r (l + a + (3; a, (3) = n(l + (3) > . (8.12) 

a— >oo 

8.3 Corollary (Bounds on the Outer Boundaries). For < (3 < 1 and 

a > 1, £/iere exzsi constants 7, which are all independent of a and (3 such that 

-7i < Mr(l + a + /9;«, / 3) < - 72 < , (8.13) 

and 

< 73 < «r(l + a + /3',a,/3) < 74 < 00 . (8.14) 

Proof. The proof for n r is essentially the same as the proof for v r , so we 
will only deal with u r . Because the limit as a — > 00 is strictly negative by the 
previous lemma there is a large a such that for a > a we have the desired 
lower and upper bounds. Next we simply use compactness for the rest of the 
strip. 
Q.E.D. 
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